﻿ 3 CONVOLUŢIA Introducem operatia de convoluţie in timp discret (suma de convolutie) si in timp continuu (produsul de convolutie) Calculul răspunsului sistemelor liniare şi invariante in timp, la un semnal de intrare oarecare http://shannon etc upt ro/teaching/ssist/Cap3 pdf Suma de convoluţie Sinteza semnalului de intrare Produsul intre un impuls Dirac intarziat cu k si semnalul x[n] extrage valoarea esantionului x[k]:   0x n n x n           x n n k x k n k               Semnalul x[n] este o suma de esantioane plasate la toate valorile k posibile:  x n x k n k      k 1 •Consideram un sistem liniar S, intrare x[n], iesire y[n] d x[n] y[n] S d Sistem liniar •Raspunsul sistemului la impulsul unitar, δ[n], sau raspunsul la impulsul unitar (functie pondere), h[n]: [n] h[n] δ S d Sistem liniar •In plus daca sistemul este si invariant in timp, raspunsul la impulsul unitar intarziat: δ[n-k] h[n-k] S d 3 SLIT Răspunsul sistemelor discrete liniare şi invariante în timp (SLITD) la un semnal de intrare oarecare                    ddy n S x n S x k n k k               dx k S n k x[n] y[n] k   h[n]  x k h n k        Sistem SLIT k   (LTI system) Suma de convolutie sau convolutia in timp discret:  y n x k h n k x n h n h n x n         k 2 Exemplul 1) Convolutia are durata finita             knykxnyx k  Doua semnale cu durata finita N si N 12 Convolutia are durata N+ N-1 12 Suportul (durata) semnalului: un semnal are N valori semnificative (posibil si unele nule), dar dupa care valorile sunt identic nule 5 Exemplul 2) Convolutia are n durata infinita x n a n               2h n n n N 1n 1ka            n N x h n a20 1 1a k   2Nn  211nNka         n N x h n a a2 1a   21k n N 6 3 Semnale digitale cauzale Daca semnalul de intrare, respectiv sistemul, sunt cauzale, semnalul de iesire este si el cauzal: n x n h n n y n x k h n k0 and 0 , 0                 0k Daca sistemul este cauzal: h n n h n h n n n Z0 , 0 ,                n y n x n k h k x k h n k                  0kk7 BIBO Stabilitatea •Sistem stabil: Dacă semnalul de intrare este mărginit atunci şi răspunsul trebuie să fie mărginit, “bounded input bounded output” (BIBO) •Conditia de BIBO stabilitate Un sistem digital SLITD este stabil daca si numai daca raspunsul sau la impuls este absolut sumabil:  1 h k ,h n l     8 k 4 Stabilitatea acumulatorului n  h n n y n x k            k   n - cauzal, , acumulatorx n y n x k        0k n  1x n n y n n n              0k Acumulatorul este instabil •Semnalul de iesire nu este marginit •Folosit in practica: Se limiteaza intervalul de insumare n , sau, din cand in cand se pune semnalul de iesire pe zero Proprietati Elementul neutru •Impulsul Dirac [n] este element neutru pentru convolutie x[n]* [n] = x[n], pentru orice semnal x[n] δ[n] h[n] h[n] h[n] – raspunsul la impuls al sistemului 10 5 Sistemul identitate •Are raspunsul la impuls δ[n] (elementul neutru pt convolutie) •Pentru intrare δ[n]  iesirea este tot δ[n] 11 Sistem de intarziere Raspunsul la impuls h[n]=[n-n] 0 Raspunsul sistemului este o varianta intarziata a semnalului de intrare      00x n n n x n n 12 6 Asociativitatea convoluţiei Conectarea în cascadă (serie) a SLIT                 nhnhnxnhnhx Convolutia este asociativa  2121          2121hhxhhxy  21hhhe Raspunsul la impuls al sistemului echivalent este       12h n h n h n           12y n x n h h n Prin conectarea în cascadă a 2 SLIT stabile se obţine tot un SLIT stabil 111             2121lnhhlnh,lnh Suma de convoluţie este comutativă  nhnhnhnhnh e1221  nx  nx  nh  nh  nh nh 1 221  ny  ny   nx  nx 1 2 La conectarea în cascadă nu contează ordinea 7 Sistemul invers Sistemul cu raspunsul la impuls este inversul sistemuluihn   i cu raspunsul la impuls hn daca prin conectarea lor in   cascada se obtine un sistem identitate h n h n n        i  Un exemplu de sistem invers Pentru sistemul de intarziere :       0nnnh Sistemul invers h n n n      0i Deplasarea in timp in sensul invers pe axa timpului Nu este un sistem cauzal Distributivitatea convoluţiei faţă de adunare Conectarea în paralel a SLITD •Convolutia este distributiva fata de adunare                    1 2 1 2x h n x h n x h h n 16 8 Demonstratie                    nhxnhxnynyny 2121                knhkxknhkx   21   kk               knhknhkx  21 k         nhhx 21 17 Functia indiciala a unui SLITD, s[n]: Răspunsul la treapta unitara n  x n n y n s n h n n h k,                     k     1h n s n s n       Daca sistemul este cauzal, n 0, pentru 0h n n s n h k           0k Functia indiciala a unui acumulator este un semnal rampa x n n s n n n    9 SLITD cu răspuns finit la impuls (FIR) şi cu răspuns infinit la impuls (IIR) FIR IIR b , 0nnM     0ahn 0 , in rest  •Raspuns la impuls cu durata finita (FIR) sau cu durata infinita (IIR) FIR •Ecuatia cu diferente finite NM     kka y n k b x n k 00kk •Un sistem FIR are proprietatea : 0 and 0a a a a     0 1 2N 20 10 Semnalul de iesire •Folosind ecuatia cu diferente finite: M kb       knxny  a k0  0 •Folosind definitia convolutiei:          knxkhny  k   21 FIR •Prin identificare se obtine: b  k  ,0M,k, 1,   kh a  0  0, restin  •Se observa ca raspunsul la impuls are durata finita, M+1, de aici si numele “sistem cu raspuns finit la impuls” 22 11 Sisteme IIR (Infinite Impulse Response) NM        aaknxbknya ; 0 si 0   kk 10 kk 00  0 5 1y n y n x n                    nhnynnx condiţie iniţială nulă 1 0;y   0 0 5 1 0 1 0 1 0 0y y x y y h                    1 0 5 0 1 0 1 0 5 1 1y y x y y h                   2 2 0 5 1 2 0 2 0 5 2 2y y x y y h                   3 3 0 5 2 3 0 3 0 5 3 3y y x y y h                   n 0 5h n n     Implementarea SLITD caracterizate prin ecuaţii cu diferenţe finite, liniare şi cu coeficienţi constanţi •SLITD sunt descrise matematic prin ecuatii cu diferente finite cu coeficienti constanti •Sunt implementate folosind subsistemele: celule de memorare (intarziere), multiplicatoare cu o constanta, sumatoare 12 Sisteme de ordinul unu, Implementarea directă I            nxbnxbnyanya 11 1010 •O celula de memorare: intrare x[n], iesire x[n-1] Doua multiplicatoare cu o constanta, b0 si b1: zn•   b0x[n] si b1x[n-1] •Un sumator Implementarea sistemului nerecursiv, de mediere alunecatoare (Moving average) MA Iesire: z[n] Implementarea sistemului recursiv, autoregresiv (Autoregressive) AR Iesire: y[n] 25            nxbnxbnyanya 11 1010 26 13 NM Sistem de ordinul N         ;  kkknxbknya kk 00  M1N y n z n a y n k                 kknxbnzk  k010ka  Celule de memorare, multiplicatoare, sumatoare 27 Forma transversala pentru FIR •Consideram un sistem FIR : 1      nzny  a 0 •Sub-sistemul 2 = multiplicator cu 1/a 0 •Substituim cele M sumatoare din subsistemul 1 cu un singur sumator cu M intrari 28 •Am presupus ca a=1 0 14 Implementarea directă II •Ordinea sistemelor nu este importanta (conexiune in serie) •Daca se schimba ordinea: se obtine o forma echivalenta implementarii directe I, si anume implementarea directa II Se inlatura sistemele de intarziere 29 redundante Sistem de ordin N, Implementarea directă II M N b, b, M+1N+2 a, a, …, a N+1N+2M b = 0 = 0 …, N 30 15 Produsul de convoluţie Răspunsul SLITC la un semnal de intrare oarecare •SLIT descris matematic de operatorul S Trebuie sa gasim iesirea y(t) atunci cand se cunoaste intrarea x(t) x(t) y(t) S SLIT (LTI system) •Reamintire: Proprietatea de filtrare a impulsului Dirac δ(t)        0d         •Daca vom considera functia test x(τ), iar ca impuls, varianta deplasata in timp cu t (aici variabila timp este τ), atunci avem:  x t d x t             31  y(t) x t x t d          S   SLIT •Semnalul de iesire      ()y t S x t S x t d x S t d                          const const      Produsul de convolutie a semnalelor x si h Se poate calcula raspunsul unui sistem cunoscut, cu raspunsul la impuls ( h ) la un semnal de intrare oarecare ( x )  y t x t h t x h t d                   32 16 Proprietati •Elementul neutru: distributia Dirac, δ(t)   x t t x t x t,          •Convolutia este comutativa aproape peste tot (a p t ) Cu notatia t – u = τ, avem:  f g t g u f t u du g f t         •Convolutia este distributiva in raport cu adunarea                    1 2 1 2x h t x h t x h h t •Convolutia este asociativa f t g t h t f t g t h t                     33 Cateva remarci •Daca x(t), h(t)  L1 , clasa functiilor absolut integrabile (sistemul este stabil), atunci y(t)=x(t)*h(t)  L1 •Daca x(t), h(t)  L2 , clasa functiilor cu patratul modulului integrabil (energie finita), atunci convolutia x(t)*h(t) exista, este marginita si continua •Daca unul din factori este din L1, iar celalalt din L2, de ex semnalul de intrare este de energie finita, x(t)  L2 si sistemul este stabil h(t)  L1, atunci iesirea este de energie finita y(t)=x(t)*h(t)  L2 34 17 Exemplul 1) semnale de durata finita          1f t t t T          2g t t t T Durata convolutiei = T+T (suma 12 duratelor celor doua semnale) 0, 0t   2,0t t T         2 2 1,f g tT T t T      1 2 1 1 2,T T t T t T T    120,t T T 35 •Se noteaza g(-τ) = z(τ) •Se deplaseaza in timp z(τ) cu t catre dreapta: 0, 0;t f g t    t       20 , ;t T f g t d t  0 t       2 1 2, ;T t T f g t d T  2tT 2T          1 1 2 1 2, ;T t T T f g t d T T t  1tT        36 12, 0, 0 t T T f f t f g t 18 Exemplul 2) TT        ;f t t t       22    ;g t f t     t   1;f g t T t T t T                T  37 Exemplul 3) 11  , ?f t f t f t         1tt •f(0); f(t) si | f(t)| sunt pare 1 •functia f apartine lui L   du          022 2arctg 2f t dt u  12u 0 •f * f este convergenta a p t; dar nu si in t=0  11          20f f f f d d  38 1    19 Exemplul 4) TT    f t t t ;           22    t  , 0 1;g t a t a       2 durata infinitaf g L   T  tT  21  1f g t a t            12    ln  a T  tT  2 1 at     2    39 • t < -T/2 , f*g(t) = 0 •–T/2 <t <T/2 suprapunere partiala: tT  1t  2t 1f g t a d a       1  lnT 2 a •t ≥ T/2 suprapunere completa: T 2TT tt   221t f g t a d a a      1 T ln40 2 a 20 Condiţia ca un SLITC să fie cauzal            de sistem cauzal 0Exemplu 0   ,tththt,th t                          dhtxdthxtxthty     0         ttxtx Daca si semnalul de intrare este cauzal, se :obtine tt                          ;dhtxdthxtxthty   00        ttyty Exemplu de sistem necauzal Conditia de BIBO stabilitate a SLITC •Stabilitate “bounded input bounded output” Un sistem continuu SLIT, este stabil daca si numai daca raspunsul sau la impuls este absolut integrabil  1 h d ,h t L           42 21 Integratorul Not bounded Răspunsul la impuls, h(t), al unui integrator este σ(t) x t t y t h t t                Răspunsul la un semnal marginit, de exemplu, σ(t), este semnalul rampa, care este nemarginit  t,t  0                         ttdtt    t, 00     Integratorul este un sistem instabil Se foloseste in practica fiindca semnalele de intrare sunt de durata finita (caz in care semnalul de iesire va fi marginit) Răspunsul indicial al unui SLITC •Raspunsul indicial=raspunsul la treapta unitara t s t h t t h d        •Derivata sa de ordinul unu este raspunsul la impuls 's t h t  •Sistem cauzal: t  s t h t t h d       0 44 22 •Raspunsul sistemului h la semnalul rampa ;x t t t y t S t t h t t                      •Derivata de ordinul doi a lui y(t) este h(t) " " y t h t t h h            x t t t      y t h t t t        S y t h t      LTI system 45 •Raspunsul la impuls h(t) este: –derivata de ordinul intai a raspunsului indicial s(t) –derivata de ordinul doi a raspunsului la semnalul rampa 46 23 Semnificatia practica a proprietatilor produsului de convolutie •Sistemul echivalent obtinut prin conectarea in serie a doua sisteme are raspunsul la impuls     2 1 1 2h t h t h t h t h t •Sistemul echivalent obtinut prin conectarea in paralel a doua sisteme are raspunsul la impuls 47       12h t h t h t Sistemul invers Sistemul identitate •Doua sisteme conectate in serie, al doilea sistem invers, atunci iesirea este semnalul de intrare original y(t)=x(t) •Sistemul echivalent, h(t) * h(t), este un sistem i identitate h(t)* h(t) = (t) i 48 24 Implementarea SLITC caracterizate de ecuaţii diferenţiale liniare, cu coeficienţi constanţi Exemple de SLITC caracterizate de ecuaţii diferenţiale liniare cu coeficienţi constanţi: derivator, integrator ()tdx (tdy ) ()ty=() tx= dtdt Forma generală a ecuaţiei diferenţiale : kk NNtxdtyd        0 Nkkkkaba  kk 00 dtdt Implementarea in forma directa II cu circuite de derivare NkNktxdtyd        0 Nkkkkaba  kk 00 dtdt •Sistem de ordinul N •Forma directa II: sistemele de intarziere (timp discret)  circuite de derivare (timp continuu) •Dificil de construit 50 25 Implementarea in forma directa II cu circuite de integrare •Sunt preferabile circuitele de integrare •Integram de N ori ecuatia diferentiala  o ecuatie integrala Cu notatiile            ,10 , ,10 , , , NNxxxyyy ,  •obtinem ecuatia integrala: NN              kNkkNktxbtya  kk 00 51 kk NN d y t d x t  0a b a k k Nkk 00kkdt dt  0,y t y t  t      111y t y t t y d    2t         1 1 22y t y t t t y d d      12kkt         1 1 2 11 kkkky t y t t y d d d d               01 ; , x t x t x t x t t  52 26 53 Exemple i) Sistemul de ordinul intai ()()NtxkdNtykd =∑ ba ∑ kkkk 0=0=kdtkdt 1=bRCaa,, = 1= 010 27 Exemplu, sistem de ordinul doi 2 d y t dy t     LC RC y t x t       2 dtdt NN ,a y t b x t      kkn k n k     00kk           0 1 2 2LCy t RCy t y t x t         •Prin identificare se obtin coeficientii, si obtinem forma II de implementare 1   0b a a RC a LC0 1 21 ; ; 55 Structura transversala pentru sisteme FIR N •Raspunsul la impuls:   kTthth k 0k  ,hNTtx hTtxhtxty N10 N  ,thtxkTthtxty k k0 N  kTthth k 56 k0 28 